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We present here a detailed multifractal scaling study for the electronic transmission resonances with 
the system size for an infinitely large one dimensional perfect and imperfect quasiperiodic system 
represented by a sequence of 5-function potentials. The electronic transmission resonances in the 
energy minibands manifest more and more fragmented nature of the transmittance with the change 
of system sizes. We claim that when a small perturbation is randomly present at a few number 
of lattice sites, the nature of electronic states will change and this can be understood by studying 
the electronic transmittance with the change of system size. We report the different critical states 
manifested in the size variation of the transmittance corresponding to the resonant energies for both 
perfect and imperfect cases through multifractal scaling study for few of these resonances. 
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I. INTRODUCTION 

The discovery of quasiperiodic phase in Al-Mn alloys 
by Sehechman et. al [Q led to a great interest in under- 
standing the generic properties of such systems from elec- 
tronic structure, transport and phonon dynamics point of 
view. The simplest possible choice being to study an one 
dimensional quasiperiodic system. Merlin and co-workers 
first studied the Fibonacci superlattices where the 
unusual fractal like structure in the spectral properties 
giving rise to a characteristic features in experiments. 
Thus the theoretical understanding of such one dimen- 
sional quasiperiodic systems and the recent advances 
of molecular-beam-epitaxy (MBE) make these systems 
ideal candidates to study and to carry out experiments 
on transport properties || of these systems. Electronic 
states of one dimensional Fibonacci quasiperiodic sys- 
tems have been attempted in the recent past | 10 1^] 
to explore some global characteristics in them. It is well 
known that for such one dimensional quasiperiodic sys- 
tems the spectrum is a cantor set and has self-similar 
structure which results from the hierarchy of Fibonacci 
lattices |i[-|TT| . Most of these previous works have mainly 
concentrated on studying critical nature of energy bands 
as well as the nature of wave function from localization 
point of view for isolated chains. But less attention has 
been given to the study of the transmittance in detail to 
look for localization and delocalization aspects through 
some rigorous scaling analysis pd| , ^2| , p2]j24| . Apart from 
this, most of the calculations have been performed start- 
ing from the tight-binding Anderson Hamiltonian where 
one models quasiperiodicity either through the diago- 
nal on-site energy or through the off-diagonal hopping 



term. However, the actual quasiperiodic potential in the 
Schrodinger's equation may not be possible to decouple 
into a diagonal term which will only retain quasiperiod- 
icity and a constant hopping term and vice versa. So 
one should expect that some of the physical features one 
finds in the numerical investigation through the conti- 
nous models may explain more closely the experimental 
situation. It is possible to fabricate semiconductor het- 
erostructure with a variety of potential profiles along the 
growth directions, such as rectangular and parabolic ones 
in resonant tunneling devices and saw-tooth potential in 
delta doped layers Apart from that, it has been re- 
alized in the recent past that in some of the thin quasi 
one dimensional wires, the variation of the potentials at 
the impurities is much larger than that of the energies in- 
volved by the weak electric field; the Kronig-Penny model 
is then the natural choice for a relevant model to start 
with. 

In recent years, many authors have attempted the 
problem of studying the scattering of free electrons by 
an array of (5-function potentials to investigate the elec- 
tronic states and transport properties of one dimensional 
disordered systems, quasiperiodic sequences, incommen- 
surate systems ]I^-|l9|] and in particular the transport 
properties of quasiperiodic superlattices Jl3| . Here we 
have attempted the study of size dependence of trans- 
mittance through the multifractal scaling analysis. We 
have focussed our attention in characterising the trans- 
mittance resonances in the different energy minibands 
that one observes in an array of one dimensional Fi- 
bonacci (^-function potentials. This has been done in 
the same spirit where in a disordered chain one stud- 
ies the transmittance vs energy plot having signature of 
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exponentially localized states along with a few stochastic 
resonances also known as Azbel resonances [ ^4|Jl5| l. 

The strong variation of transmittance with system size 
gives an indication of the non-trivial nature of trans- 
mittance for the various resonant energies. Based on a 
detailed numerical investigation of the fluctuating pat- 
tern of the transmittance within the framework of mul- 
tifractal scaling of transmittance with system size, it is 
possible to infer on the derealization and localization 
aspects for both perfect and imperfect one dimensional 
quasiperiodic systems respectively. In this work we will 
see that both the perfect and imperfect realisations re- 
semble some novel size variation which is intermediate 
between the so called exponential decay and the homo- 
geneous Bloch like contributions. We claim that many 
of these are resonant states are of nearly critical nature 
instead of showing either exponential decay or homoge- 
neous spread along the chain. We have made a rigor- 
ous study for the electronic transmittance in imperfect 
Fibonacci system to verify whether a minor change in 
the potential in a perfect chain can destroy the critical 
nature of the states altogether leading to a complete lo- 
calization of electronic states. Recently, some works in 
the area of imperfect quasiperiodic chains have been re- 
ported [^3| claiming the instability of the spectrum of a 
Fibonacci superlattice due to introduction of a single de- 
fect. We have seen that for many of the localized states 
the strong spatial fluctuations in the transmittance will 
exhibit nearly self-similar features like critical states upto 
a large length scale followed by a long-tailed oscillatory 
pattern. However, many of the resonant critical states 
due to the presence of imperfectness also change but show 
some fluctuations reflecting some kind of instability in the 
a-f(a) spectrum subject to large size variation. 



II. BRIEF OUTLINE OF THE MODEL 

Let us consider the dynamics of electrons in a Fibon- 
naci quasiperiodic array of (^-function potentials. We 
know that a Fibonacci sequence Sj of order j is obtained 
by j successive applications of the transformations : 



Sj+i 
A 
B 



Sj © Sj. 

AB 

A 



(1) 



Where Sq = A and S\ — AB and the symbol © means 
the operation of composition. For a very large but finite 
j one can generate an infinitely large array of 5-function 
potentials in the Fibonnaci sequence where the following 
relations hold: 



lim 



F j+1 = Fj + Fj-i and 
Fj-i\ fy/%-1 



F; 



(2) 
(3) 



Fj being the number of atoms in the sequence Sj. Re- 
cently Sanchez et. al. Q have made an explicit derivation 
of the generalized Poincare map and calculated the elec- 
tronic transmittance due to scattering from a lattice with 
Fibonacci [^) and correlated disorder potentials within 
the Kronig-Penny model || . In this paper we have based 
our numerical computations of electronic transmittance 
using the same Poincare map analysis having a very sim- 
ple and straight forward mathematical structure. This 
algorithm has been numerically verified to be very stable 
and suitable for doing calculation for very large system 
size in one dimensional systems. We illustrate here only 
a brief outline of the model. A detailed derivation of the 
generalized Poincare map etc. will be found in the ref. 
||. We concentrate on the Fibonacci model for an elec- 
tron interacting with lattice potentials of the following 
form. 



N 



v{x) = 2J r„ 6{x - x n ) 



(4) 



We have chosen the position of the (5-function potentials 
to be regularly placed (i.e x n = na,a = 1) with the 
strength Ta and Tb for A and B type of atom respec- 
tively. The Schrodinger's equation is then given by ( 
using h 2 = 2m =1), 



JV 



-v 2 + ^2 r„ 5(x - X n ) 



n=l 



(5) 



One can now introduce reflection and transmission am- 
plitude Rn,Tn respectively through the following rela- 
tionship: 



V(x) = e ikx + R N e~ ikx 
$(x) = T N e lkx 



x < 1 
x > N 



(6) 
(7) 



The tranmission amplitude T/v after N sites is given by 
the recurrence relation, 

A ! i a *N-l@N y Pn \ A ,„S 

A N = (aw H 7, )An-x - (-5 )A N -2 (8) 
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(9) 
(10) 

(11) 



The initial conditions Aq = 1 and A\ = a,\ arc suffi- 
cient to determine the transmission amplitude Tn com- 
pletely for an arbitrary size i.e for a of large number of 
5-function potentials. The transmittance t/v is given by 



tjv — T n Tn — — — — — 



1 
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III. MULTIFRACTAL SCALING ANALYSIS FOR 
THE VARIATION OF TRANSMITTANCE WITH 
SYSTEM SIZE 

In the recent past, Thakur et. al. made some electron 
localization studies in an aperiodic potential where the 
variation of transmittance with the system size has been 
investigated p2|. It has been established that this vari- 
ation of the tranamittance with the increase of system 
size in an infinitely large system exhibits some intrinsic 
global features which can be fully captured by the multi- 
fractal scaling analysis ^(J. The same scaling aspects 
have also been used for the Azbel resonance problem 
[pTstl in one dimensional random chain as well as in the 
study of extended and the critical states in the vicinity 
of dimer resonance energy in the random dimer problem 
p4j . Following this line of approach, our motivation in 
the present work is to examine the electronic transmis- 
sion of different resonant characters seen in the transmit- 
tance versus energy diagram of a quasiperiodic Fibonacci 
superlattices. We observe that it is possible to explain 
the dependence of the transmittance with the system size 
through the multifractal scaling aspects of the transmit- 
tance. Here we outline the mathematical content of the 
multifractal scaling for electronic transmittance only be- 
cause of its direct connection with our computation and 
its relation to some of the typical transport features man- 
ifested through some global properties in the different 
energy minibands in our calculations. We will discuss 
the same in the context of the multifractal scaling al- 
gorithm developed by Chabbra and Jensen in brief 
This algorithm has been successfully used in analyzing 
the crossover states in generalized Aubry model, random 
dimer model and also in many other contexts in the lit- 
erature p2|j2^j25| ]. The multifractal formalism has been 
developed essentially to describe the statistical proper- 
ties of some measure in terms of its distribution of the 
singularity spectrum /(a) corresponding to its singular- 
ity strength a. In our calculation we take the normalized 
transmittance Pi as the required measure which is given 
by 



Pi = 



Ti 



l~,i=l 



where Ti is the transmittance from one end of the sys- 
tem upto the i-th segment when the length of the chain is 
divided into N equal small segments such that the trans- 
mittance for a given size is obtained by always increasing 
the previous size by adding the same number of atoms. 
According to Chabbra and Jensen if we define the q-th 
moment of the probability measure Pi by the following 
expression: 

Mi(<7, AO = 



l^i=i r i 

then a complete characterization of the fractal singulari- 
ties can be made in terms of the function fii(q, N). Hence 



one can derive the spectrum of fractal singularities de- 
fined by, (a, /(a)). The expression for f(a) is given by, 



1 ^ 



N — >oo logiV ^— ' 

6 i=l 



m(q,N) log M i(<7, N) (13) 



and the corresponding singularity strength a is obtained 

as 



a = lim 



1 



JV 

\ " 



N^oo logiV 4^ 



\it log Pi 



(14) 



In this regard, the specific nature of electronic trans- 
mittance, reflected through a max i.e large negative q, 
oi m in *- e for large positive q, and also through f(a max ) 
and f(a m i n ) values, can be inferred by the following char- 
acteristics for large N(N — * oo) : 



1. Extended nature : a m i n 

Q-max * 1; f i^^max) * 1- 

2. Localized nature : a m i n 

^rnax f (^max ) 1- 



1, f{otmin) -* 1, 
0, f(a min ) — » 0, 



3. Critical nature : a versus f(a) curve will show a 
tendency to converge onto a single curve as N — > 
oo. 

Ideally for extended nature the a versus /(a) will 
contract eventually to (1, 1) whereas for critical na- 
ture the a versus /(a) curve, with the increase of 
system size, may sometimes show minor fluctua- 
tions as well. 



IV. RESULTS AND DISCUSSION 

In this section, we discuss the results for the electronic 
transmittance of both the perfect and the imperfect Fi- 
bonacci superlattices through solving the Eq. (|5|) making 
use of the mathematical relations Eqs.(p[jl^). Here we 
mainly focus our attention on the electronic transmis- 
sion resonances in different energy minibands and have 
investigated whether the resonances survive when we in- 
crease the system size drastically by plotting the trans- 
mittance versus energy plots for various system sizes. As 
stated earlier, our aim is to examine whether the delo- 
calized characteristics of the resonances show some in- 
teresting universal properties for such quasiperiodic po- 
tentials even when subject to a large size variation. We 
finally claim that our numerical studies of multifractal 
scaling for the variation of transmittance with system size 
capture most of the generic features of electronic states 
and transport in the Fibonacci quasiperiodic superlattice 
within the Kronig-Penny model. 



3 



A. Results for a perfect Fibonacci sequence 

In this subsection we discuss the results for a per- 
fect Fibonacci quasipcriodic potential. In Fig.l we have 
shown the transmittance T versus E plots for = 2 
units and Tb = I units for three different system sizes, 
N = 75025 (lower curve), N = 196418 (middle curve), 
and N — 3 178 11 (topmost curve). One can clearly see 
that as the system size increases transmission envelope 
fragments more and more. It's also true that at any stage, 
in general, the T versus E plot has a fragmented shape. 
These characteristic features are due to the critical na- 
ture of the energy spectrum of the Fibonacci quasiperi- 
odic system. We now choose arbitrarily a transmission 
resonance corresponding to a small region of energy mini- 
bands. This resonance can be identified more clearly if 
the T versus E plot within this regime can be checked 
for a very fine energy mesh ~ 1CP 7 or even less such that 
a resonanance within this fine energy resolution will ex- 
hibit it's signature in the T versus N plot. In Fig. 2(a) 
we have chosen E = 20.0918599999603 units as one of 
such resonant energies and have shown the variation of 
the transmittance T with system size N . The T versus 
N plot clearly shows highly oscillatory pattern having a 
very complicated structure and even for variation of size 
from 10 5 to 7 x 10 5 no. of atoms or layers it retains the 
same global shape in it's pattern. We have carried out 
the multifractal scaling of the transmittance with it's sys- 
tem size to analyze this spatial pattern in Fig. 2(c). Also 
in Fig. 2(b) we have shown another plot of T versus N for 
2 x 10 5 to 3 x 10 5 no. of atoms or layers for a resonant 
energy E = 3.204406000 units. We see that the oscilla- 
tory pattern in Fig. 2(b) seems to be different from that of 
Fig. 2(a) altogether. In Fig. 2(c) the a — f(a) plot shows 
a small inward contraction with the increase of system 
size from N = 196418 to N = 317811 no. of atoms or 
layers and then from N = 317811 to N = 514229 no. of 
atoms or layers. The deviation of a m in and a max values 
with the size variation is much less and is <~ 0.01. We 
identify them by saying nearly critical resonant states. 
In Fig. 2(d) the a versus f(a) plot corresponding to the 
variation of the transmittanceT with system size N in 
Fig. 2(b), have been shown. Here we see that as the sys- 
tem size has been increased from N — 75025 no. of atoms 
or layers to N = 196418 no. of atoms or layers the curve 
deviates from the initial position as a whole but when 
the system size has been increased to N = 317811 no. of 
atoms or layers, the curve comes in between the two pre- 
vious positions. However, here the amount of deviation 
of a max values are ~ 0.025 to 0.040 and a m i n ~ 0.01 to 
0.02. One knows that for a localized state with the in- 
crease of systen size a max values systematically goes out 
but here for the size N = 317811 it comes inside. This 
is the signature of critical state showing that it's nei- 
ther localized nor extended but have self-similar pattern 
for the size variation upto an infinitely large system size. 
We identify this resonant state as critical resonant state 



. We now consider set of parameters for the potentials, 
namely, = 3.0 units and Tb = 1.0 units and check 
whether here also one finds some typical resonant states. 
In Fig. 3(a), the T versus N plot has been shown for 
the resonant energy E = 6.07038530468009 units. The 
variation of the transmittance (T) with system size (N) 
manifests a highly fragmented and oscillatory behavior 
from 5 x 10 5 to 832040 no. of atoms or layers. We also 
choose another resonant energy E — 15.0532885305118 
units for the same potential parameters Ta and Tb where 
the T versus N plot shows clearly some oscillatory as 
well as fragmented pattern. However, here, the overall 
shape approxomately gives a periodic oscillation. The 
f(a) versus a plot corresponding to the variation of 
transmittance(T) with system size(A) in Fig. 3(a) has 
been shown in Fig. 3(c). Note that as we increase the sys- 
tem size from N = 75025 to N = 317811 both the a max 
and a min values change by very small amount ~ 0.01 or 
less. As we increase the size to N — 514229 , the curve 
begins to deviate outwards but when we increase the size 
from N = 514229 to N = 832040 no. of atoms or layers 
the curve then again shows an inward contraction. So 
we see here that the qualitative behavior is much like a 
resonant state in Fig. 2(c) and this can be identified as 
a critical resonant state as well. Next we analyze the a 
versus f(a) plot in Fig. 3(d) for the resonant state corre- 
sponding to the energy E = 15.0532885305118 units for 
which the variation of transmittance with system size has 
been shown in Fig. 3(b). Here we notice that, as we in- 
crease the system size from N = 196418 to N = 317811 
no. of atoms or layers the curve undergoes an inward 
contraction. Next with the subsequent increase of sys- 
tem sizes upto N = 1346269 no. of atoms or layers, the 
curves systematically contracts as we increase the system 
size more and more, the curves show a relative tendency 
to contract by very lesser amount. We identify this to be 
qualitatively similar to the resonant states as depicted in 
Fig. 2(c). So we again identify it to be a nearly critical 
resonant state and it is in some sense different from the 
critical resonant shown in Fig. 2(b) and Fig. 3(c). 

We think these two typical resonant states, namely, 
critical and nearly critical can be identified for different 
choice of potentials also. It's also true that they are of 
two distinct nature which one can differentiate through 
the a versus /(a) plots. These resonant states exist for 
different energy regimes and also for any choice of 
and Tb- This particular aspect is true since we have 
chosen the resonant energies in Fig. 2(a) and Fig. 2(b) or 
in Fig. 3(a) and Fig. 3(b) in completely different energy 
regimes arbitrarily. 

B. Results for an imperfect Fibonacci sequence 

Now we go over to analyze the effects of disorder 
present at random in any layer. The amount of pertur- 
bation in the potential parameters Ta and T b have been 
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considered 5% at random in any layers through perturb- 
ing 2% of the atoms or layers in the system as a whole. 

In Fig. 4(a) we have shown the transmittance(T) ver- 
sus system size(iV) plot for the resonant energy E = 
6.07062311828030 units corresponding to starting Fi- 
bonacci potentials Ta — 3.0 units and Tb = 1-0 unit 
and then a fractional change in the potntials /= 0.05 
have been considered randomly at any layer where 2% 
of the layers have been perturbed as a whole. One can 
see that the transmittance amplitudes look like different 
patches of spikes of many different heights. The variation 
of transmittance with system size(iV) has been shown 
from 10 6 to 1346269 no. of layers. Note that it looks nei- 
ther like extended nature or exponentially localized but 
appranctly the profile of the spikes seem to have some 
power law decay. We have checked through the multi- 
fractal scaling analysis whether the resonant states would 
show some systematic decay or not. 

In Fig. 4(b), the plot of f(a) versus a corresponding 
to the variation of transmittance with system size in 
Fig. 4(a), have been shown. We have shown a versus /(a) 
curves corresponding to the size variation starting with 
the initial size 832040 no. of atoms or layers and then 
succesively increasing no. of atoms or layers by 10 5 no. 
of atoms or layers. One can clearly see in Fig. 4(b) that 
a versus f(a) curves show some instability even for the 
variation of transmittance beyond 10 6 no. of atoms or 
layers. This instability in the a versus f(a) plot is due 
to the effect of disorder in the Fibonacci quasiperiodic 
layers. However, if one looks at the overall variation of 
the transmittance with system size, from 10 6 to 1346269 
no. atoms or layers then it shows the decay but of vary 
complicated form and this is also captured in the outer- 
most a versus f(a) plot corresponding to the variation 
of transmittance T for the system of size from 832040 to 
1346269 no. of atoms or layers eventually. One can see 
that this is the typical long-tailed fluctuations which we 
have observed due to the interplay of the effect of disor- 
der and the resonant nature of critical state in the ideal 
system. Since this is a state which shows a size variation 
in between critical and exponentially localized state, one 
can identify this resonant state, due to the presence of 
imperfectness in the system, to be more closely like crit- 
ical resonant state. Next we have shown, in Fig. 5(a), the 
variation of transmittance with system size for a resonant 
state corresponding to the energy E = 4.3987600006124 
units and for the values Ta — 2.0 units, Tb = 1.0 units. 
It's a typical localized state with fluctuations throughout 
but having less fluctuations in the tail part as compared 
to that in the beginning, i.e, around 50, 000 to 2 x 10 5 no. 
of atoms or layers. In Fig. 5(b), the a versus f(a) plots 
corresponding to the variation of transmittance with sys- 
tem size as shown in Fig.5(a), for N = 317811,514229 
and, N = 832040 no. of atoms or layers, have been 
shown. It's quite clear that with the increase of system 
size the curves begin to move outwards reflecting the very 
nature of a localized state having an overall systematic 
decay. In Fig. 5(c) we have shown the plot of transmit- 



tance with the system size N for an imperfect superlattice 
with E = 20.88552401746 units. This is also a resonant 
state with strong fluctuations throughout which arises 
because of the presence of disorder in the system. In 
Fig. 5(d) we have shown the a versus /(a) plot corre- 
sponding to the variation of the transmittance(T) with 
system size (iV)in Fig. 5(c). Note that, as the system size 
has been increased from N = 196418 to 317811 no. of 
atoms or layers, the curves undergo an outward shift but 
as we increase the system size from 317811 to 514229 
no. of atoms or layers it again goes in between the two 
previous positions. This sort of instability reflected in 
the a versus f{a) is due to the presence of imperfect- 
ness coupled with the long-ranged quasiperiodic order in 
the system. Since these states do not show any system- 
atic decay either in the a — f(a) plot, and show some 
kind of critical nature, we identify them more closely like 
critical resonant states also. Therefore, to paraphrase, 
the manifestation of the effect of disorder on the reso- 
nant critical states, apart from the existence of having 
systematic decay, has been found to be twofold: firstly, 
it has long-tailed fluctuations and secondly having sus- 
tained strong fluctuations even when subject to large size 
variations as mentioned above. 



V. CONCLUSION 

In this work we explore the nature of electronic states 
and transport in one dimensional ^-functional Fibonacci 
potentials from a rigorous study of electronic transmis- 
sion resonances corresponding to the different energy 
minibands. In this regard, although some studies of 
transmittance have been reported in the past, our nu- 
merical study of the size dependence of the same is di- 
rectly associated with the characterization of the nature 
of electronic states and transmittance using the multi- 
fractal scaling aspects which to our knowledge has not yet 
been investigated adequately till date. We should point 
out the fact that many of the resonances show strong os- 
cillations with the variation of the system size and the 
oscillating pattern continues to survive till the width of 
the minibands, in which it belongs, goes to zero and an 
energy gap opens up there. In a perfect Fibonacci case, 
critical behavior of the transmittance with system size 
has its origin from the typical nature of the Fibonacci 
potential and at least a few of them can be distinguished 
from the others. More specifically, here, the clear distinc- 
tion between the different critical resonant states have 
been possible with the help of multifractal scaling analy- 
sis. On the other hand, we have considered the imperfect 
situation to be a case where the imperfectness appears 
randomly as a homogeneous substitutional disorder in 
the different lattice points causing only minor changes in 
the Fibonacci potentials which may be present e.g in an 
experimental realisation of Fibonacci superlattice. 

In the imperfect case we have examined that this 
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may induce localization effects which may be manifested 
through some non-exponential spatial decay of transmit- 
tance and the resonant states in the perfect situation will 
appear again more closely like critical resonant states but 
with some modified nature. 

The long-tailed oscillatory pattern of the transmit- 
tancc shows a very slow spatial decay and it retains its 
nearly self-similar pattern even if one drastically changes 
the system sizes. Although the overall pattern of the os- 
cillations eventually show decaying fluctuations, the de- 
cay is not uniform— indicating an interplay between the 
long-ranged quasipcriodic order and the disorder an elec- 
tron experiences with the increase of system size. The 
other manifestation of the interplay between long-ranged 
quasiperiodic order and substitutional disorder is the in- 
stability of the resonant state in the variation of trans- 
mittance with the increase of system size. This has been 
also reflected in the a-f (a) spectrum in our numerical 
studies. We, therefore, hope that it will be interesting 
to search for both non-exponential localization and crit- 
ical properties of the electronic transmittance through 
designing some controlled experiments in this direction. 
We believe that our work will shed some more light in 
understanding the theoretical notion of the novel scaling 
properties for the electron localization and derealization 
aspects directly in the size dependence of the conduc- 
tance in Fibonacci quasiperiodic systems. 
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FIG. 2. (a) Transmittance T versus system size N for a 
pure Fibonacci lattice with Ta ~ 2 units, Ts = 1 unit and 
E = 20.0918599999603 units. 

(b) Transmittance T versus system size N for a pure Fi- 
bonacci lattice with Ya = 2 units, Yb = 1 unit and 
E = 3.20440600 units. 

(c) f(a) versus a plot for the same set of parameters as in 
Fig. 2(a) but for different system sizes having maximum no. 
of atoms or layers N = 196418, (line), N = 317811 (dashed 
line) and for N = 514229 (smaller dashed line), where in each 
case starting size is 5 no. of atoms or layers and the rest is 
obtained increasing 5 no. of atoms or layers each time. 

(d) f(a) versus a plot for the same set of parameters as in 
Fig. 2(b) but for different system sizes having maximum no. 
of atoms or layers N = 75025, (line), N = 196418 (dashed 
line) and for N = 317811 (smaller dashed line), where in each 
case starting size is 5 no. of atoms or layers and the rest is 
obtained increasing 5 no. of atoms or layers each time. 

FIG. 3. (a) Transmittance T versus system size N for a 
pure Fibonacci lattice with Ta = 3 units , Y B = 1 unit and 
E = 6.07038530468009 units. 

(b) Transmittance T versus system size N for a pure Fi- 
bonacci lattice for the same set of potential parameters as 
in Fig.3(a) but for energy value E = 15.0532885305118 units. 

(c) f(a) versus a plot for the same set of parameters as in 
Fig. 3(a) for different system sizes having maximum no. of 
atoms or layers N = 75025 (line), N = 317811 (dashed line) 
, N = 514229 (smaller dashed line) and N = 832040(dotted 
line), where in each case starting size is 5 no. of atoms or 
layers and the rest is obtained increasing 5 no. of atoms or 
layers each time. 

(d) f(a) versus a plot for the same set of parameters as in 
Fig. 3(b) for different system sizes having maximum no. of 
atoms or layers N = 196418 (line), N = 317811(dashed line) 
, N = 514229 (smaller dashed line), N = 832040(dotted line) 
and for N — 1346269(dash-dotted line), where in each case 
starting size is 5 no. of atoms or layers and the rest is obtained 
increasing 5 no. of atoms or layers each time. 



FIG. 4. (a) Transmittance T versus system size N for an 
imperfect Fibonacci lattice for Fa = 3.0 units ,Tb = 1 unit for 
E = 6.07062311828030 units. The fractional change / = 0.05 
has been considered randomly at any layer or atom (say i) 
so that the new r, = Yi + f * Yi. The concentration of such 
defects present in the lattice is c = 0.02 

(b) /(«) versus a plot for the same set of parameters as in 
Fig. 4(a) for different system sizes. The data set have been 
considered starting with size having no. of atoms or layers 
N = 832040 and adding 10 5 no. of atoms or layers (line), 
adding 2 x 10 5 no. of atoms or layers (dashed line), 3 x 10 5 no. 
of atoms or layers (smaller dashed line), 4 x 10 J no. of atoms 
or layers (dotted line) and 5 x 10 5 no. of layers (dash-dotted 
line) with steps of 5 layers in each case. 

FIG. 5. (a) Transmittance T versus system size N for an 
imperfect Fibonacci lattice for the set of parameters Fa = 2.0 
units and Y B = 1.0 unit for E = 4.39876000006124 units . The 
fractional change / = 0.05 has been considered randomly at 
any layer or atom so that the new T; = r» + / * Yi . The con- 
centration of such defects present in the lattice is c = 0.02. 

(b) /(«) versus a plot for the same set of parameters as in 
Fig. 5(a) for different system sizes. The data set have been 
considered starting with size having maximum no. of atoms 
or layers N = 317811 (line), N = 514229 (dashed line.) and 
N = 832040(smaller dashed line) no. of atoms or layers, where 
in each case starting size is 5 no. of atoms or layers and the 
rest is obtained increasing 5 no. of atoms or layers each time. 

(c) Transmittance T versus system size iV for an imperfect 
Fibonacci lattice for the same set of potential parameters as 
in Fig.5(a) but for E = 20.88552401746 units. 

(d) f(a) versus a plot for the same set of parameters as in 
Fig. 5(c) for different system sizes. The data set have been 
considered starting with size having maximum no. of atoms 
or layers N = 196418 (line), N = 317811 (dashed line) and 
N = 514229 (smaller dashed line.), where in each case start- 
ing size is 5 no. of atoms or layers and the rest is obtained 
increasing 5 no. of atoms or layers each time. 

[26] 
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